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SUMMARY 

Water  losses  from  canals,  pools,  etc.  , exposed  to  the  effects  of 

nearby  nuclear  explosions  are  caused  by  ground  motions,  alrblast  loading, 

* 

and  impacting  debris.  Certain  of  these  phenomena  are  reasonably  well 
understood  and  can  be  modeled  with  confidence,  e.g.,  fluid  motion  Induced 
by  ground  motion  and  overpressure  loading.  Other  phenomena  that  could 
contribute  to  fluid  loss  are  less  well  understood  and  can  be  modeled  only 
with  poor  confidence,  e.g.,  the  growth  of  surface  instabilities  under 
dynamic  pressure  loading.  Impacting  debris  is  representative  of  a third 
class  of  loss  mechanisms  in  that  the  water  loss  could  conceivably  be 
prescribed  as  a function  of  Impact  velocity,  mass  and  shape  of  debris, 
etc.;  yet  still  be  poorly  characterized  because  of  the  uncertainties  in 
specifying  the  debris  environment.  This  report  omits  any  analysis  of 
water  loss  produced  by  impacting  debris  (although  debris  certainly  con- 
tributed to  the  fluid  loss  observed  on  high-explosive  cratering  experi- 
ments) . Rather,  we  concentrate  on  fluid  motion  Induced  by  ground  motion 
and  overpressure  loading — and  provide  some  comment  on  fluid  response 
associated  with  complex  interactions  between  the  alrblast  dynamic 
pressure  and  the  fluid's  free  surface. 

Our  approach  is  to  assume  the  existence  of  a causal  relationship 
between  the  maximum  wave  response  and  the  water  loss,  and  to  obtain  data 
to  quantify  this  relationship.  Analytic  solutions  for  fluid  response 
(in  rectangular  containers)  to  ground  motion  excitation  are  obtained  and 
compared  to  data  from  controlled  laboratory  shake-table  experiments.  We 
find  that,  for  horizontal,  sinusoidal  ground  motion  loading  conditions 
of  frequency  f,  the  depth  removed  (Ad)  is  proportional  to  the  amplitude 
of  the  displacement  input  (a),  i.e.,  Ad  ■ \(f)*a  where  \(f)  is  a fre- 
quency-dependent constant  of  proportionality.  The  theoretically  deter- 
mined maximum  wave  response  (h)  is  also  proportional  to  the  amplitude 
of  the  input  displacement  amplitude;  i.e.,  h = Y(f)*a  where  Y(f)  is  a 
frequency-dependent  constant  of  proportionality.  Furthermore,  it  appears 
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that  Y(f)  ^ rx(f)  where  T Is  a constant  for  all  rectangular  geometries 
examined.  Based  on  this  observation,  we  hypothesize  that  the  fractional 
depth  loss  Is  given  by 

M - rJl 
d ■ ‘d  ’ 

where  F Is  a constant  for  a given  container  geometry  (F  0.31  for  rec- 
tangular containers)  and  d is  the  container  depth.  The  above  hypothesis 
is  used  to  estimate  fluid  loss  from  containers  exposed  to  combined  alr- 
blast-horlzontal  ground  motion  environments.  Such  predictions  are  con- 
sistent with  data  obtained  from  MIXED  COMPANY,  although  the  estimates 
lie  on  the  low  side  of  the  data  scatter.  This  disagreement  may  be  due 
to  the  effects  of  dynamic  pressure,  vertical  ground  motion,  and  impacting 
debris,  which  are  ignored  in  the  wave  height  calculation. 

Following  the  encouraging  comparison  with  MIXED  COMPANY  data,  we 
estimate  fluid  losses  from  typical  large-scale  containers  subjected  to 
alrblast  and  ground  motions  predicted  for  nuclear  explosions.  This 
analysis  suggests  that  Che  relative  portion  of  fluid  loss  associated 

with  ground  motion  varies  considerably  with  yield,  container  dimension 

1" 

an.  geologic  variations  of  interest.  In  some  cases,  airblast  loading 
dominates.  For  example,  when  large  trenches  (^300-ft  length)  in  dry 
soil  are  subjected  to  low-yield  surface  bursts,  airblast  effects  dominate 
the  close-in  fluid  loss,  while  ground  motion  effects  dominate  the  close- 
in  fluid  loss  for  trenches  exposed  to  large  yields  in  wet  soil.  The 
worst  case  is  associated  with  wet  geologies 

These  results  may  underestimate  the  actual  fluid  loss  because  con- 
tributions associated  with  dynamic  pressure  and  debris  impact  have  not 
been  considered.  Additional  studies  are  required  to  assess  the  water 
loss  likely  to  be  connected  with  such  effects.  However,  it  is  believed 
that  only  minor  perturbations  to  the  fluid  loss  predictions  presented 
here  will  be  introduced  by  these  phenomena. 


’^'Alrblast'induced  water  loss  is  independent  of  geology  except  in  a 
second  order  sense  which  considers  coupled  airblast  and  ground  motion 
effects . 
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Water  loss  from  trenches,  canals,  pools,  etc.  , may  be  prevented  or 
significantly  reduced  with  the  use  of  berms  placed  around  their  perimeters. 
As  a practical  matter,  the  cost  of  berms  would  be  nominal  because 
excavated  soil  would  be  readily  available  for  their  construction.  The 
waveheight  calculations  discussed  in  this  paper  could  be  used  to  develop 
design  criteria  for  berm  features  (e.g.,  berm  height)  once  the  threat 
environment  and  soil  geology  are  specified. 
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SECTION  1 
INTRODUCTION 


The  effects  of  nuclear  bursts  on  water-filled  canals,  pools,  etc., 
are  of  Interest  for  a variety  of  reasons  [1].  The  analysis  presented 
here  is  concerned  with  airblast  and  ground  shock  induced  fluid  motion 
in  a rectangular  trench,  where  the  blast  wave  is  normal  to  the  trench 
sidewalls  (see  Figure  1.1).  We  emphasize,  however,  that  the  choice  cf 
this  particularly  simple  geometry  is  motivated  only  by  our  desire  to 
simplify  notation.  Most  of  the  results  obtained  here  can  be  easily 
generalized  to  more  complicated  situations — both  in  the  sense  of  con- 
tainer geometry  and  in  the  sense  of  blast  wave  direction  relative  to  the 
container  orientation.  In  any  case,  the  geometry  Illustrated  in  Figure  1.1 
is  an  excellent  approximation  to  nearly  all  the  field  and  laboratory 
experiments  discussed  In  Section  5.  Since  the  information  obtained 
from  these  experiments  consists  almost  entirely  of  post-shock  fluid  loss 
mersurements , the  theoretical  effort  presented  here  is  directed  primarily 
toward  the  development  of  a suitable  fluid-loss  prediction  technique, 
although  information  on  other  important  fluid  response  characteristics 
(e.g.,  Che  velocity  field)  is  also  obtained. 

In  principle,  a complete  description  of  fluid  response  in  containers 
exposed  to  blast  and  shock  effects  from  nuclear  explosions  can  be  derived 
from  the  governing  equations  for  the  appropriate  boundary-initial  value 
problem.  Because  ground  displacement  waveforms  are  mainly  composed  of 
low  frequency  components,  which  are  essentially  unaffected  by  the  presence 
of  shallow-burled  structures,  the  free-fleld  ground  motions  can  be  used 
as  a displacement  boundary  condition  for  the  pool  side  walls.  In  addi- 
tion, the  effect  of  airblast  overpressure  can  be  accurately  represented 
as  a pressure  boundary  condition  for  the  fluid's  free  surface.  However, 
although  the  ground  motion  and  overpressure  boundary  conditions  can  be 
specified  precisely,  interactions  between  the  airflow  (dynamic  pressure) 
and  the  free  surface  are  not  well  understood  and  can  only  be  approximately 
modeled.  This  fact,  along  with  our  suspicion  that  the  water  loss  caused 
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by  dynamic  pressure  related  forces  may  be  of  minor  importance  for  the 
larger  length  scales  of  interest,  lead  us  to  first  consider  the  effect 
of  ground  motion  and  overpressure  on  fluid  response  in  Section  2,  where 
a discussion  of  the  complete  boundary-initial  value  problem  leads  to 
scaling  rules  employed  to  design  HE  field  experiments;  and  in  Section  3, 
where  a detailed,  linearized  analysis  of  fluid  motion  is  presented.  An 
examination  of  dynamic  pressure  effects,  including  scaling  requirements, 
is  presented  in  Section  4,  while  in  Section  5,  we  consider  the  available 
experimental  data  in  formulating  a fluid  loss  prediction  technique. 
Finally,  Section  6 provides  estimates  of  fluid  loss  in  large-scale  con- 
tainers subjected  to  alrblast  and  horizontal  ground  motion  loading. 

The  water  loss  induced  by  Impacting  debris  can  also  be  modeled 
theoretically  provided  the  debris  size,  shape  and  impact  velocities  are 
known.  We  shall  not  attempt  to  treat  the  effect  of  debris  impact  here. 
Other  studies,  beyond  the  scope  of  this  paper,  will  consider  possible 
bounds  to  water  loss  associated  with  impacting  debris  (as  well  as  the 
potential  hazards  of  Impacting  missiles  on  the  survivability  of  the 
contents  of  liquid  shelters) . 
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SECTION  2 

GOVERNING  EQUATIONS  AND  SCALING  FOR 
GROUND  MOTION  AND  OVERPRESSURE -EFFECTS 

^ ^ Scaling  Requirements 

An  examination  of  the  governing  equations  for  fluid  response  in 
open  containers  indicates  that  at  least  five  separate  scaling  require- 
ments must  be  satisfied  to  guarantee  dynamic  similarity  of  the  fluid 
motion,  In  a practical  sense,  these  requirements,  which  are  associated 
with  gravitational,  viscous,  surface  tension,  cavitation  and  compressible 
effects,  cannot  be  simultaneously  Incorporated  into  the  design  of  model 
experiments.  It  appears,  however,  that  elimination  of  the  non-gravita- 
tlonal  scaling  requirements  may  not  seriously  affect  the  properties  of 
fluid  motion  which  are  important  to  blast  and  shock-induced  fluid  loss. 

In  particular,  viscosity  and  surface  tension  probably  have  a very  minor 
effect  on  parameters  such  as  maximum  wavehelght  and  net  fluid  loss,  pro- 
vided container  dimensions  are  greater  than  a few  Inches.  Likewise, 
cavitation  (which  may  not  even  occur)  and  fluid  compressibility  are  only 
important  during  passage  of  a strong  air  shock  (>100  psi)  over  the  fluid's 
surface,  Air  shock  related  compressive  effects  are  currently  being 
investigated;  however,  for  the  purposes  of  this  report  we  shall  assume 
that  they  do  not  seriously  affect  the  fluid  response  characteristics 
mentioned  above. 

Consider  an  open,  two-dimensional,  rectangular  fluid-container 
system^  of  length  Z and  mean  fluid  depth  d.  Assume  that  the  free  sur- 
face of  the  fluid  is  subjected  to  an  overpressure  loading  P(x,t)  and 
that  the  container  undergoes  a rigid  body  displacement,  defined  by  the 
motion  of  surrounding  ground  material,  D(t).  By  neglecting  the  fluid's 


T 

The  specialization  to  two-dimensional,  rectangular  containers  is 
made  only  to  simplify  the  notation.  Generalization  of  the  following 
procedure  to  a three-dimensional  system  of  arbitrary  geometry  is  straight- 
forward. The  results  are  unchanged. 
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viscosity,  surface  tension,  and  compressibility,  and  by  introducing  the 

length  scale  Z and  time  scale ■\^£/ g,  the  equations  of  motion  and  boundary 

1* 

conditions  can  be  expressed  in  the  non-dimensional  form 

2 ^ 

= 0 (Continuity  equation),  " ^ ^ 


8n 

3t 


* 

+ v*n  • 


* 

it 

3z  * * 

2 =n  . 


<2. 2) 


tt 


P(x,t) 
* Pg^ 
n 


* 


+ n 


= 0, 


(2.3) 


.tt 


* 

3(t> 

* 

3x 


x*=$(t)  ‘t/«. 


9x 


x*=l+b(t)»i/S, 


1 dP(t)  . 

I dT  * ^ 


(2.  A) 


tt 


3 

i 

dZ 


1 dP(t) 
~ dx 


“>■ 

k 


z =-(d+D(t)-k)/«. 


where  gravity  acts  in  the  negative  z-direction,  (f  is  the  velocity  potential 

n is  the  displacement  of  the  free  surface  from  equilibrium,  1 and  k are 

* * * 

unit  vectors  in  the  x and  z-directlons  and  x =x/A,  z =zlZ,  n =r\/Z, 

For  a discussion  of  the  dimensional  equations,  see  Reference  2. 
tt 

These  equations  represent  the  kinematic  free  surface  boundary  con- 
dition (2,2),  the  dynamic  free  surface  boundary  condition  (2,3),  and  the 
solid  wall  boundary  conditions  (2,4). 
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T=t/V”i-777 


* * * ■ *• 

dx  dz 


To  guarantee  similitude  in  Equations  2.1  through  2.4,  we  must  require 


P(x,t)  = PgUp  (x  ,t) 


(2.5) 


D(t)  =2,D  (t). 


(2.6) 


Thus,  if  a fluid-container  system  of  characteristic  length  , is  sub- 
jected to  an  overpressure  Pj^(x,t)  and  a displacement  Dj^(t),  a geometri- 
cally similar  fluid-container  system  of  characteristic  length  will 
undergo  dynamically  similar  motion  provided 


D2(t)  = 


p2(x,t)  = 


(2.8) 


m 


With  viscous,  surface  tension,  and  compressible  effects  neglected,  Equa- 
tions 2.7  and  2.8  completely  specify  the  ground  motion  and  overpressure 
scaling  requirements,  and  in  the  idealistic  situation  where  complete 
control  over  the  environment  is  feasible,  no  more  need  be  said.  Unfor- 
tunately, complete  environmental  control  is  not  always  possible,  especi- 
ally when  HE  field  tests  are  used  to  model  prototype  response.  For  example, 
although  HE  and  NE  ground  motion  waveforms  are  reasonably  similar  (and 
hence  Equation  2.7  can  be  satisfied  by  adjusting  the  length  scale  and  range 
from  ground  zero),  the  spatial  and  temporal  dependence  of  overpressures 
generated  by  explosions  of  different  yields  is  inconsistent  with  the 


t 

The  position  x = 0 is  defined  as  the  leading  edge  of  the  container. 
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If  the  overpressure  positive  phase  duration  is  small  compared  to 
the  period  of  the  lowest  fluid  response  mode  and  effects  of  the  negative 
phase  are  negligible,  the  analysis  presented  in  Section  3 (see  Equation  3.13) 
demonstrates  that  with  regard  to  airblast-induced  fluid  motion,  positive 
phase  impulse  gradient  is  the  most  Important  overpressure  parameter. 
Integrating  Equation  2.8  over  time,  we  obtain  the  scaling  relation 

I^Cx)  = . (2.9) 

where  I denotes  positive  phase  overpressure  Impulse.  If  P is  the  peak 
overpressure  at  x = 0,  R is  the  radius  from  ground  zero  and  W is  the 
explosive  yield;  impulse  scaling  is  obtained  through  the  following  steps: 

1/2  1/3 

Since  I 'Vi  P W provided  P ^ 1000  psi  (see  Ref  3),  Equation 
2.9  requires 

The  geometrical  requirement  specified  by  the  function  arguments 
in  Equation  2.9  is  satisfied  if  R 'v-  Jl. 

The  relation  P (W^^^/R)'^  {where  2 ^ n ^ 3 for  P ^ 10  psi) 
together  with  (1)  and  (2)  implies  that  P 'v  Z W™,  where 
j = l/(l+2/n)  and  m = 2/3(3-J),  i.e.,  1/2  5.  j 1 3/5  and 
5/15  < 5/18. 

1/3 

For  fixed  P,  P W and  the  impulse  gradient  ("vI/P)  exhibit  a negli- 
gible yield  dependence.  Thus,  the  only  explicit  scaling  rule  is  P 'v  , 
although  there  is  an  implicit  yield  dependence  because  we  have  assumed 
that  the  positive  phase  duration  is  small  compared  to  the  lowest  mode 
period.  For  very  large  yields  and/or  very  small  containers,  this  type 

of  scaling  is  inappropriate  (e.g.,  impulse  scaling  is  accurate  if 
1/3 

W /T<<1,  where  W is  expressed  in  MT  and  T is  the  lowest  mode  period  in 
seconds).  Since  the  fundamental  periods  for  containers  with  lengths  of 


1. 

2. 

3. 


100-300  ft  vary  from  5-15  sec,  the  use  of  impulse  scaling  for  yields 

5<>1  MT  could  lead  to  Inaccuracies.  For  situations  where  this  condition  is 

1/4 

violated,  and  the  yield  is  not  sufficiently  large  for  W scaling  (see 
below)  to  apply,  or  when  the  effects  of  a negative  phase  are  of  Interest, 
no  simple  scaling  rule  applies,  and  a more  detailed  examination  of  the 
problem  is  necessary. 

2 . 3 Quarter  Root  of  Yield  Scaling 

1/4 

The  scaling  requirements,  P "v  A and  I V , commonly  referred  to  as 

quarter  root  of  yield  scaling,  are  derived  in  Reference  4 with  the  aid 

3 1/2 

of  the  strong  shock  relations  W — PR  and  u(shock  speed)  ~P  , under 

the  assumption  that  the  characteristic  time  scale  for  the  problem  is  the 

air  shock  transit  time  across  the  container. 

If  the  positive  phase  duration  is  much  longer  than  the  container's 
1/4 

lowest  mode  period,  W scaling  is  applicable  because  the  important 
forces  are  applied  to  the  fluid  surface  only  during  the  air  shock  tra- 
verse (and  hence  the  transit  time  does  serve  as  a characteristic  time 
scale).  On  the  other  hand,  in  cases  where  the  positive  phase  duration 
is  short,  important  forces  are  applied  to  the  surface  long  after  the 

air  shock  front  has  passed  (when  the  overpressure  decay  time  constant 

1/4 

is  unrelated  to  the  shock  transit  time).  Unfortunately,  W scaling 
appears  to  be  appropriate  only  for  extremely  large  yields  and/or  very 
small  containers.  For  container  lengths  on  the  order  of  a few  hundred 
feet,  this  scaling  is  not  accurate  unless  the  yield  K,  100  MT. 

2 . 4 Ground  Motion  Scaling 

For  a given  D^^(t),  Equation  2.7  defines  a 0^(1) , which  will  preserve 
dynamic  similarity  of  the  fluid  motion.  Laboratory  experiments  can,  in 
principle,  satisfy  Equation  2,7  to  any  desired  degree  of  accuracy;  therefore, 
prototype  waveforms  can  be  accurately  modeled — to  the  extent  that  they 
are  known.  From  the  standpoint  of  field  tests,  ground  motion  scaling 
is  approximate  because  of  difficulties  in  precisely  scaling  HE  and  NE 
ground  motion  waveforms.  However,  a well-designed  HE  experiment  provides  a 


wide  variety  of  ground  motion  environments  which  usually  model  some 
nuclear  waveforms  of  interest. 


If  the  dominant  frequency  of  a prototype  displacement  waveform  is 

f , and  the  dominant  frequency  of  a field  test  waveform  is  f , the 
P t 

length  scale  which  should  be  used  In  the  field  test  is  determined  from 

Equation  2.7  as 


i ~ U I ’ 

p \ p; 


where  la  the  prototype  length  scale. 

Thus,  for  field  test  experiments,  ground  shock  scaling  requirements 
can  be  satisfied  fairly  accurately,  but  only  at  the  expense  of  the  over- 
pressure requirements.  This  condition,  coupled  with  the  fact  that  rig- 
orous scaling  of  overpressure  effects  is  often  impossible,  strongly  sug- 
gests that  we  attempt  to  understand  fluid  response  from  first  principles, 
rather  than  totally  relying  on  information  from  approximately  scaled 
experiments.  Consequently,  in  the  next  section,  we  present  a more 
detailed  analysis  of  fluid  response  to  ground  shock  and  overpressure. 
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SECTION  3 

DYNAMICS  OF  FLUID  RESPONSE  TO  GROUND  MOTION  AND  OVERPRESSURE 
3.1  Formulation 

In  this  section,  a normal  mode  analysis  is  used  to  examine  the 
problem  illustrated  In  Figure  1.1.  Since  the  technique  oE  modal  decom- 
position can  be  applied  to  any  fluid-container  system  excited  in  an 
arbitrary  direction,  the  extension  of  these  rectangular  trench  results 
to  more  complicated  situations  is  straightforward. 

Consider  a two-dimensional,  rectangular  container  of  length  i 
in  the  x-direction  in  a gravitational  field  g which  acts  in  the  negative 
z-directlon.  Suppose  that  the  equilibrium  fluid  depth  in  the  container 
Is  d and  that  the  sides  of  the  container  are  extended  sufficiently  to 
prevent  any  fluid  loss . 

With  viscosity,  surface  tension,  and  compressibility  neglected. 
Equations  2.1  through  2.4  and  the  initial  conditions  b(o)=n(o)=0  comprise 
a boundary-initial  value  problem  with  a (presumably)  unique  solution. 
Unfortunately,  the  highly  nonlinear  nature  of  the  problem  precludes  the 
use  of  simple  solution  techniques,  although  it  is  possible  to  obtain  a 
useful  approximate  solution  by  linearizing  Equations  2.2  through  2.4 — i.e., 
we  assume  that  the  free  surface  and  ground  displacements  are  small  compared 
to  container  dimensions  and  that  the  particle  velocities  are  small  compared 
to  typical  surface  wave  speeds.  These  assumptions  are  justified  for  most 
cases  of  interest,  and  even  in  situations  where  this  is  not  so,  the  line- 
arized solution  is  at  least  qualitatively  useful. 

The  linearized,  dimensional  versions  of  Equations  2.2  through  2.4  are 


(x.t) 


(3.1) 


^ (x,o,t)  + g77(x,t)  + 


(3.2) 


II  (o.z.t)  - II  a,z,t) 


(x,-d,t)  = H(t) 


(3.3) 


which  together  with  ^^(11=0  and  the  aforementioned  initial  conditions  com- 
prise the  linearized  problem.  Here,  D(t)  and  H(t)  respectively  are  the 
horizontal  and  vertical  components  of  the  ground  displacement. 

In  terms  of  defined  by 

<|i  = i)(t)(x-|-)  + H(t)z  + ^ 


the  governing  equations  become 

= 0 


(x,o,t)  + (g+H(t)>»j(x 
ot 


,t)  + + D(t)(x-  -|)  = 0 


' I?  (x.t) 

II  (o.z.t)  = II  (i.z.t)  = ol 


(3. A) 


(3.5) 


(3.6) 


(3.7) 


II  (x.-d.t)  - 0 J • 

Since  the  crater-induced  vertical  accelerations  are  usually  small  compared 
to  g,  ii(t)  will  be  omitted  from  Equation  3.5  for  the  remainder  of  this 
section.  Section  4 discusses  the  effect  of  this  term  under  more  general 
circumstances.  Note  that  H(t)  in  Equation  3.6  only  requires  the  equilib- 
rium position  of  the  free  surface  to  move  vertically  with  the  displacement 
H(t). 

Equations  3.4  through  3.7.  with  all  vertical  motion  terms  omitted 
can  be  satisfied  by  setting 
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•n 


B (t)cos 
n 


where 


A (t)  + gB^(t)  * D(t)«  + jS„(t) 

X\  II  11  11 

- B„(t) 


with  a and  6 defined  by 
n n 


(x-i/2) 


oo 

- / n cos  (n7rx/j2  ) 


n=l 


(i,e.,  a^  = 41/ (im)^  for  odd  n;  - 0 for  even  n) • 

OO 

P(x,t)  = )3^(t)COE(n77X/j2) 

n=l 

(l.e.,  e^(t)  = ^ ^ P(x.t)cos  dx). 

Elimination  of  A from  Equations  3.8  and  3.9  yields 
n 


2 

B + u»  B = 
n n n g 


K + %] 


where 

1 /2 

(j  = [(n7rg/j2)tanh(nTrd/j2)] 


Is  the  natural  radian  frequency  of  normal  mode.  Equation  3.10  can 

be  Interpreted  as  describing  the  displacement  of  a harmonic  oBcillator 

2 

of  frequency  u , which  is  excited  by  the  force  (w  /g) (Da  +6  ).  The 
’ n n n n 

solution  of  Equation  3.10,  which  satisfies  the  initial  conditions 

B (o)=B  (o)"0,  is 
n n 


B„(t) 


2 

a a 

n n 


D(t)  - 0) 


D(t)  sin  (hi  t 
n 


t-T])dTi 


f 

Jo 


(3.11) 


)sln(o)^tt-T])dT  , 


where  t“0  marks  the  start  of  the  initial  disturbance  (either  ground  motion 
or  overpressure).  By  differentiating  Equation  3.11  and  substituting  the 
results  into  Equation  3.9,  one  can  calculate  A^(t)  and  thereby  determine 
the  fluid  velocity  field,  Vifi 

The  extension  of  Equations  3.10  and  3.11  to  fluid-container  systems 
with  more  complicated  geometries  is  accomplished  by  replacing  with 
the  actual  natural  frequencies  of  the  system  in  question  and  the  modal 
functions  cos(nirx/Jl)  by  the  true  mode  shapes,  which  may  be  functions  of 
two  horizontal  space  variables.  Note  that  this  replacement  changes  the 
definitions  of  a^  and  In  laboratory  experiments,  the  first  few 

can  be  measured  easily,  whereas  the  mode  shapes,  with  the  possible  excep- 
tion of  the  first,  are  probably  very  difficult  to  determine.  Fortunately, 
giound  shock  simulation  experiments  (to  be  discussed  in  Section  5)  indi- 
cate the  existence  of  a good  correlation  between  fluid  loss  and  the  maxi- 

•f* 

mum  virtual  waveheight  associated  with  the  first  mode  alone.  Thus, 


Since  the  calculation  of  B (t)  assumes  that  fluid  cannot  leave 
the  container,  these  wavehelghts* are  termed  virtual. 
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for  fluid  loss  prediction,  the  first  ntode  response  appears  to  be  the 
roost  Important. 

If  ground  motion  effects  are  ignored.  Equation  3.11  reduces  to 


B„(t) 


r)sln(w^(t-r))drdx h <3.12) 


3 . 2 Fas t Decay  Approximation 


If  the  modal  periods  are  long  compared  to  the  positive  phase  duration 
(l.e. , (2Tt/a)^)>>D'*^) , Equation  3.12  can  be  approxloated  by 


-2ai  slnCcLi  t)  -i  . V 

= —ja • 


where  I(x)  is  the  positive  phase  Impulse,  and  the  effects  of  a negative 
phase  have  been  Ignored  (the  negative  phase  Impulse  is  much  smaller  than 
the  positive  phase  impulse  when  P ^ 100  psi).  Since 

/ \3/2  - 


= 1 - 


I(R+x)  . 1 

^P(R+x)l 

I(R)  ' 

^ P(R)  J 

is  a good  approximation  for  yields  and  overpressure  of  interest  [2], 
Equation  3.12  can  be  further  approximated  by 

r / /I  \ 2 


■ sinCw^t),  n odd 


B (t)  = 
n 


0,  n even. 


(3.13) 


Equation  3.13  indicates  that  the  surface  response  is  dominated  by  over- 
pressure impulse  gradient  (~I(R)/R).  If  impulse  gradient  scaling  rules 
are  employed,  dynamic  similarity  of  the  fluid  motion  Is  guaranteed  pro- 
vided (27r/u  )>>D^. 
n 

3.3  Step  Wave  Approximation 

In  the  other  extreme,  2ir/u)^<<D^,  and  the  overpressure  can  be  approxi- 
mated by  a step  function: 
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I; V 


P(x.t)  = 


O 0<t<x/u 


P t>x/u  , 

O 


(3.14) 


where  u is  the  air  shock  propagation  velocity.  Substitution  of  Equation 
3.14  into  Equation  3,12  gives 

f P - t /k  \ 


" (4)  # « 


B^(0  = 


(3.15) 


0,  n even 


where  T = Jl/u  is  the  shock  transit  time  and  T “ 2tt/o)  . The  similitude 
t 1/4  " 

properties  of  Equation  3.15  are  consistent  with  W scaling,  which  requires 

1/2  1/2 

P 'V  s,  and  assumes  that  u "v  P £ . 
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Sections  2 and  3 have  been  concerned  with  the  somewhat  idealized 


problem  of  overpressure  and  ground  motion  induced  fluid  response.  Clearly, 
the  high-speed,  time-dependent  air  flowfield  associated  with  airblast 
environments  can  also  Interact  with  the  fluid  free  surface,  exciting 
waves  and  possibly  causing  fluid  lo.  « in  the  process.  Although  this 
interaction  is  difficult  to  describe  theoretically,  progress  can  be 
made  in  assessing  the  significance  of  some  of  the  more  simple  dynamic 
pressure  effects. 

A qualitative  understanding  of  dynamic  pressure  effects  can  be 
obtained  by  ignoring  many  of  the  more  complex  phenomena  which  would 
otherwise  render  the  problem  intractlble.  For  example,  if  we  ignore 
the  boundary  layer  and  the  Inherent  three-dimensionality  of  the  flow- 
field,  treat  the  problem  as  quasi-static,  and  assume  that  the  deformation 
of  the  fluid  surface  is  small,  the  resultant  airflow-fluid  interaction 
can  be  interpreted  as  a Kelvin-Helmholtz  instability.  These  assumptions 
generally  exaggerate  the  strength  of  the  airflow-fluid  interaction.  In 
particular,  boundary  layer  effects  tend  to  reduce  the  effective  air 
velocity;  large  free  surface  deformations  induce  flow  separation  and 
thereby  limit  the  growth  of  surface  instabilities;  and  the  actual 
three-dimensional  flow  can  relieve  itself,  laterally,  around  the  wave 
disturbance . 

Consider  a steady,  inviscid  airflow  over  a small  amplitude,  sinusoidal 
surface.  For  subsonic  conditions,  the  surface  variations  induce  a pres- 
sure disturbance  in  the  air  given  by 


8p  - 


(«.l) 


where  q is  the  amplitude  and  k the  wave  number  of  the  surface  deflection, 

Q is  the  unperturbed  dynamic  pressure  associated  with  the  airflow,  and 
s 
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Is  a drag  coefficient  that  depends  on  the  flow's  Mach  number  [5]. 
Following  an  approach  used  to  describe  airblast  interactions  with 
structures  [6],  we  assume  that  Equation  4.1  can  be  extended  to  unsteady 
flow  situations  by  setting 

SP(x.t)  = -C^(M^)Q^(t)r;(x,t)k  (4.2) 

where  Qg(t)  is  now  the  time-dependent  dynamic  pressure,  and  is  an 

average  drag  coefficient,  assumed  to  depend  only  on  the  initial  (t  = 0) 

Mach  number  M . 

o 

The  pressure  fluctuations  (4.2)  can  be  Incorporated  into  (3.5)  by 

replacing  P(x,t)  with  P(x,t)  + 6P(x,t).  The  net  result  of  this  procedure 
2 

is  to  replace  m B in  Equation  3.10  with 
n n 


d s 
1^8 


B 

n . 


Thus,  as  expected,  the  dynamic  pressure  effect  (represented  by  Kelvin- 
Helmholtz  instability)  forces  the  modal  frequencies  to  become  time- 
dependent  and  possibly  imaginary.  As  noted  previously,  vertical  ground 
accelerations  cause  a similar  effect. 

4.2  Approximate  Treatment  of  Kelvln-Helmholtz  Wave  Generation  and 
Vertical  Ground  Motion 

The  discussions  presented  In  Section  4.1  and  3.2  indicate  that 
both  Kelvin-Helmholtz  instability  and  vertical  ground  motions  modify 
the  basic  modal  equations  as  follows: 


B (t)  + u 
n n 


(1  + ^ [t])B  (t) 
n n 


2 r 


D(t)  a + 
n 


6^(t)| 


2 


g 


(4.3) 


where 


C ,Q  (t)  / \ 

^ 1^1  (Kelvin-Helniholtz) 


e^(t)  =• 


Denoting  as  the  solution 


(Vertical  Ground  Motion) , 


2 

<»Jr.  r- 


subject  to  B (o)=B  (o)=0,  we  note  that  the  differential  equation  (4.3), 
n n 

together  with  these  same  initial  conditions,  is  equivalent  to  the 
Volterra  integral  equation 


B (t)  = B 
n 


(T)sin(w  [t-T]c’T. 
n 


The  Neumann  series  solution  to  Equation  4,4 


B (t)  = B 
n : 


(T)sin(w  [t-TjdT 
n 


+ cu  I C (T)sln(D 
n I n I 


fr  (t')B*(t')sin(u  [T-t'])dt'( 
n I 'n  ^ ^ 


L 1.  ' 


Is  useful  If  the  dynamic  pressure  and  vertical  motion  perturbations  are 
small,  In  which  case 

r'  * 

B^(t)  = B^(T)^^(T)sln(u^[t-x] )dx,  (4.6) 


where  B (t)  Is  given  by  Equation  3,11, 
n 

4.3  Application  to  the  Kelvln-Helmholtz  Mechanism 

Usually  Q (t)  can  be  approximated  by  [3] 

. . . s 


Q^<t)  = 


In  many  cases  of  Interest,  0 is  much  smaller  than  the  periods  of  the  first 
few  normal  modes,  so  that  for  the  purposes  of  integration 


sin(oJ  [t-T])2:aj  (t-r) 
n n ’ 


and  B (r)  can  be  approximated  as 
n 


B (T)  = TV  , 
n’ 


where  is  a characteristic  velocity.  Equation  4.6  can  then  be  approxi- 
mated by 


roo 


B (t)  = V t + V ^d  ^ 
n n n 


(t-T)Te”^'^®dr 


(4.7) 


provided  CO ^t«l  and  t»9.  The  relative  importance  of  the  Kelvln-Helmholtz 
correction  can  be  expressed  in  terms  of  the  parameter 


(4.8) 
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Thus,  if  K <<1,  the  series  solution  Equation  4.5  is  useful,  and  the  influence 
n 

of  the  K-H  mechanism  is  minor;  if  K »1,  Equation  4.5  probably  does  not 

n 


converge,  and  the  K-H  effect  is  probably  significant.  For  most  prototype 


environments  K <<1,  at  least  for  the  first  few  modes,  which  are  most 
n 


important  with  regard  to  fluid  loss. 


Since  (n9/i)  , K-H  wave  generation  becomes  increasingly  more 


Important  for  the  higher  modes,  (until  limited  by  surface  tension  or 
viscosity)  and  large  6 (l.e.,  large  yield)  environments.  Relatively 
speaking,  the  Importance  diminishes  as  container  dimensions  Increase. 
This  discussion  of  K-H  wave  generation  ignores  many  important  phenomena 
including  boundary  layer  effects.  Miles  [7]  has  suggested  a method  of 
evaluating  the  K-H  effect  for  steady,  incompressible,  boundary  layer 
flow.  As  one  might  expect  intuitively,  his  results  Indicate  that  the 


determination  of  should  not  involve  the  free  stream  velocity,  but 


rather  a velocity  more  typical  of  the  boundary  layer.  This  effect  would 
tend  to  reduce  the  interaction  strength  about  one  or  two  orders  of 
magnitude.  Miles'  theory  explains  in  part  why  the  K-H  mechanism  is 
unimportant  for  low-speed,  wind-water  Interactions, 

Assuming  one  wished  to  study  the  K-H  mechanism  experimentally, 


proper  scaling  can  be  obtained  by  holding  the  parameter  K^  constant . 


This  can  be  accomplished  by  requiring 


Q ,M  = constant 
s o 


(4.9) 


A, 


Since  0 - W for  explosive  tests,  the  second  relation  in  Equation  4.9 
implies  that 


Before  leaving  the  subject  of  dynamic  pressure  effects,  we  should 
point  out  that  fluid  can  be  removed  from  pools  without  the  appearance 
of  large  amplitude  waves.  Shock  tube  experiments  [8]  have  indicated 


m 


i 


i' 
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that  the  air  flow  behind  a shock  wave  can  significantly  erode  the 
fluid  surface  via  a shear  Interaction  between  the  air  and  fluid. 
It  appears  that,  under  more  realistic  field  test  explosive  condi- 
tions, shear  erosion  is  much  less  significant.  In  any  case,  a 
suitable  physical  model  of  this  effect  has  not  yet  been  perfected 
4,4  Application  to  Vertical  Ground  Motion 

Substitution  of  = H(T)/g  into  Equation  4.5  gives 


B (t)  B (t) 
n n 


- ^ B* 

nj  g n 


(t)  sin(o)  [t-T])dT. 


Often  H(t)  and  D(t)  are  non-zero  only  in  a time  interval  which  is  short 
compared  to  2T\/m  , at  least  for  the  first  few  modes.  Under  such  circum 
stances,  a good  approximation  to  is  (see  Equation  3.11) 


B (t)  = (w  a/g)D(t). 
n n 


So  that  Equation  4,10  can  be  approximated  by 


2 

0)  Ct 

B (t)  ^ ^ D(t) 
n g 


I 

>1  J 8 


(t)  (t-T)dT  [ , 


Note  that  the  relative  contribution  to  B (t)  from  vertical  motion 

n 

depends  on  the  phasing  between  the  horizontal  and  vertical  components. 
The  relative  importance  of  vertical  motions  is  roughly  represented  by 
the  size  of  the  parameter 


G . f . 


where  T is  the  duration  of  the  ground  motion  and  H is  some  "average" 
vertical  acceleration.  Since  w T<1,  the  magnitude  of  typical  vertical 


msam 


accelera'tlons  must  be  somewhat  larger  than  g if  the  vertical  motion 
contribution  to  is  to  be  significant. 

A. 5 Summary 

This  analysis  suggests  that  both  vertical  ground  motion  and  K-H  wave 
generation  may  be  unimportant,  particularly  for  large  length  scales. 
Estimates  of  K (see  Equation  4.8)  for  typical  prototype  parameters  and  low 
mode  numbers  indicate  that  K^<<1,  especially  if  Miles'  correction  is 
applied.  Unfortunately,  the  K-H  mechanism  is  only  a first-order  approxi- 
mation to  the  actual  airflow-fluid  interaction.  Although  there  is 
reason  to  believe  that  under  the  highly  transient  conditions  generated 
by  a nuclear  burst,  K-H  coupling  probably  overestimates  the  strength  of 
the  interaction,  more  effort  is  needed  to  clarify  this  point.  Likewise, 
more  effort  should  be  devoted  to  the  shear  erosion  effect  described  at 
the  end  of  Subsection  4.3,  and  other  wave  generation  mechanisms  which  are 
beyond  the  scope  of  this  paper. 

An  examination  of  Equation  4.11  for  typical  ground  motion  waveforms 
indicates  that  the  vertical  motion  correction  is  usually  small.  This 
observation  appears  to  be  consistent  with  the  results  of  ground  shock 
simulation  experiments,  which  are  discussed  in  Section  5. 


SECTION  5 

FIELD  TEST  AND  LABORATORY  GROUND  SHOCK  SIMULATION  EXPERIMENTS 
5.1  Introduction 

Recently,  several  field  and  laboratory  experiments  have  been  performed 
to  study  blast  and  shock-induced  fluid  motion  in  small-scale  containers. 

All  containers  were  completely  full  at  the  start  of  each  test,  and  the 
water  loss  was  measured  after  each  test. 

On  MIXED  COMPANY  III,  a 500-ton  tangent  HE  test,  small-scale,  fluid- 
filled  containers  (hereafter  referred  to  as  the  MIXED  COMPANY  containers) 
were  exposed  to  alrblast,  ground  motion,  and  debris  environments  of  varying 
severity  [8].  These  containers,  which  are  illustrated  in  Figure  5.1,  were 
placed  ac  distances  from  ground  zero  ranging  from  300  to  900  ft  (corres- 
ponding to  250-10  psi  peak  overpressures) . The  results  and  preliminary 
Interpretation  of  these  experiments  are  given  in  Reference  9,  and  therefore 
will  not  be  discussed  here  except  to  compare  the  results  with  fluid  loss 
predictions. 

Fluid  loss  observed  on  MIXED  COMPANY  III  was  induced  by  the  combined 
action  of  airblast,  ground  shock  and  debris  Impact.  To  aid  in  unfolding 
the  relative  influenwe  of  these  water  loss  mechanisms,  laboratory  experi- 
ments have  been  conducted  to  study  fluid  response  to  ground  motions  without 
the  interference  of  other  possible  loss  mechanisms.  These  experiments 
were  conducted  on  a duax-axis  hydraulic  dynamic  simulator,  at  Wyle 
Laboratories'  Norco,  California  facility  [10],  The  primary  test  series 
exposed  water-filled  MIXED  COMPANY  containers  to  uni-axial  (horizontal) 
and  simultaneous  bl-axial  (horizontal  and  vertical)  excitation.  A test 
series  employing  somewhat  larger  cont''iners  with  several  different  wall 
slopes  was  also  conducted.  Emphasis  was  placed  on  parametric  studies 
that  varied  input  waveform  frequency,  amplitude,  duration  and  vertical- 
horizontal  phasiiig.  The  experiments  involving  combined  vertical-horizontal 
loading  were  somewha..  complicated  by  the  fact  that,  for  moderate  vertical 
motions,  the  water  separated  from  the  container  and  entered  a gravitational 
free-fall  stage.  Under  a combined  airblast  ground-shock  environment,  fluid 
separation  probably  would  be  prevented  or  at  least  inhibited  by  the  forces 


associated  with  the  alrblast  overpressure.  It  Is  Important  to  note  that 
the  theoretical  approach  to  vertical  motion  effects  presented  In  Subsection 
4.5  cannot  accommodate  fluid-container  separation.  In  any  case,  our 
efforts  to  develop  a fluid  loss  prediction  technique  have  been  guided 
mainly  by  the  single  (horizontal)  axis  test  results;  consequently,  only 
those  results  will  be  examined  In  detail. 

5 . 2 Interpretation  of  Laboratory  Experiments 

From  Equation  3.11,  with  n - 1,  6^  = 0,  the  maximum  "virtual"  wave 
length  h associated  with  the  lowest  mode  Is 


h = — tanh 

TT 


MAX 

'0<t<T 


joCt)  - 2TTf^  J D 


(t)  sin  (w^tt“T])dT/  (5.1) 


where  £ is  the  container  length  in  the  direction  of  excitation,  D(t)  is 
the  input  horizontal  displacement  waveform,  T is  the  excitation  duration, 
and  is  the  temporal  frequency  of  the  container's  lowest  normal  mode. 

The  Wyle  experiments  employed  a displacement  waveform  approximately 
represented  by 


a sin  (2TTft) 


0<t<k/f 


t>k/f 


where  k = 1,  2 or  3.  For  the  waveform  expressed  in  Equation  5.2,  Equation 
5.1  can  be  expressed  as 


h = a - tanh  Fj^(f^/f) 


(5.3) 


where  F,  must  be  determined  numerically.  A plot  of  F, (f  /f)  and  F«(f  /f) 

K.  1 o Z o 

is  presented  in  Figure  5.2.  Note  that  the  differences  between  the  two 
functions  are  negligible  except  near  resonance.  The  large  values  of  F2 
calculated  near  resonance  should  be  interpreted  carefully  because  non- 
linear, wave-growth-limiting  mechanisms  were  not  considered  in  the  analysis 
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presented  in  Section  3.  Consequently,  although  the  two-cycle  fluid  loss  is 
expected  to  be  larger  than  the  one-cycle  loss  when  difference 

Is  expected  to  be  substantially  less  than  Indicated  by  Figure  5.2.  This 
subject  Is  somewhat  academic  for  large  amplitude  ground  motions  near  ground 
zero  because  they  tend  to  have  fewer  cycles — often  less  than  a complete 
cycle. 

We  expect  that  the  observed  depth  loss  Ad  can  be  related  In  some 
manner  to  the  maximum  wave  height  achieved  In  response  to  the  displacement 
excitation.  Thus,  we  shall  correlate  the  fluid  loss  data  with  the  virtual 
wave  height  in  Equation  5.3.  For  convenience,  both  the  depth  loss  and 
wave  height  are  normalized  by  the  mean  depth  d.  Figure  5.3  plots  experi- 
mentally determined  values  of  Ad/d  as  a function  of  the  theoretical 
virtual  wave  height  h/d  (l.e.,  ta/d)[4/iT]  tanh  [itd/l]  Fj^[f^/f])  for  a 
single  cycle  displacement  pulse  (Equation  5.2  with  k ■>  1) . These  data  are 
from  tests  involving  four  container  geometries  (broadside  and  end-on 
MIXED  COMPANY  configurations,  and  broadside  and  end-on  configurations  for 
an  8"  X 15"  x 2"  rectangular  container)  with  a/d  = 0.1,  0.2,  0.3,  0,4  for 
the  MIXED  COMPANY  containers,  and  a/d  = 0.5  for  the  other  containers.  In 
addition,  several  values  of  the  input  frequency  f were  Investigated:  f = 1, 
1.5,  2.5  and  10  Hz  for  the  MIXED  COMPANY  containers,  and  f ■ 1,  2,  and  5 Hz 
for  the  8"  x 15"  x 2"  containers.  The  lowest  mode  natural  frequencies  for 
these  containers  are  2.37  and  1.55  Hz  (broadside  and  end-on,  MIXED  COMPANY) 
and  1.55  and  0.90  Hz  (broadside  and  end-on,  8"  x 15"  x 2"  container). 

To  within  the  experimental  scatter.  Ad/d  appears  to  be  proportional 

to  h/d.  This  linear  dependence  is  even  more  apparent  if  the  data  points 

where  f /f>l  are  excluded,  as  shown  in  Figure  5.4.  For  container  lengths 
o 

on  the  order  of  100  ft  or  larger  and  yields  ^ 50  Ml,  f^/f*^!.  Thus,  the 
linear  relationship  between  Ad/d  and  h/d  appears  to  be  quite  accurate  in 
the  nondlmenslonal  frequency  domain  of  most  practical  concern. 

On  the  basis  of  this  data,  we  hypothesize  that  the  fractional  depth 
loss  is  related  to  the  maximum  wave  height  via 

Ad/d  = rh/d, 

where  T =0.31+0.02  for  rectangular  containers  and  f > f , 


Figure  5.2.  Maximum  Wave  Height  Form  Factor  as  a Function  of 
Nondimenslonal  Displacement  Waveform  Frequency 
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To  Illustrate  the  Implications  of  Equation  S.4,  a more  detailed  view 
of  the  laboratory  water  loss  data  from  rhe  MIXED  COMPANY  containers  (for 
the  displacement  pulse  In  Equation  5.2,  with  k ■■  1)  is  given  In  figures  5.5 
and  5.6.  Also  shown  In  these  figures  are  predictions  based  In  Equation  5.4: 
with  h given  by  Equation  5.3  and  T ■ 0.31. 

The  disagreement  near  resonance  probably  results  from  two  effects. 

First,  nonlinear  effects  limit  wave  height  growth  and  tend  to  flatten  the 

peak  for  relatively  large  values  of  a/d.  Second,  as  fluid  is  lost  from  a 

container,  the  resonance  frequency  f lowers— an  effect  that  tends  to 

o 

broaden  rhe  resonance  peak  over  a wider  range  of  frequencies. 

Experiments  were  also  conducted  with  8"  x 15"  x 2"  rectangular  con- 
tainers that  could  be  modified  by  Inserts  to  vary  the  wall  slopes.  Although 
too  few  data  were  obtained  to  confidently  determine  the  relation  between 
water  loss  and  input  amplitude  as  a function  of  frequency  and  wall  slope, 
the  results  appear  to  be  compatible  with  the  fluid  loss  prediction  formula 
(Equation  5.4)  where  T decreases  monotonlcally  with  decreasing  wall  slope. 

The  natural  frequency  of  the  first  mode  also  appears  to  decrease  monotoni- 
cally.  However,  the  variation  In  T and  f^  appears  to  be  rather  minor  for 
wall  slopes  greater  than  18  degrees.  For  the  minimum  slope  tested  (18*), 
f^  was  lowered  by  about  15  percent  from  the  natural  frequency  for  rectangular 
containers,  while  T decreased  at  most  by  25  percent  from  0.31  to  about  0.25, 

5 , 3 Comparison  with  field  Test  Data 

The  maximum  "virtual"  wave  height  associated  with  the  lowest  normal 
mode  can,  in  principle,  be  calculated  for  any  combined  ground  motion  air- 
blast  environment.  Equation  5.4,  which  was  derived  from  experiments 
involving  only  ground  motion  loading,  may  then  be  assumed  to  predict 
fluid  loss  for  coqibined  alrblast  ground-shock  environments.  Although 
this  assumption  has  not  been  verified  through  a comprehensive  experimental 
program,  it  is  intuitively  appealing  and  will  be  used  In  the  remaining 
analysis  in  this  paper.  Beciuse  dynamic  pressure  and  vertical  ground 
motion  add  complexity  to  any  maximum  wave  height  calculation,  and 
because  these  effects  are  thought  to  be  relatively  unimportant,  the 
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Figure  5.4  Fractional  Depth  Loss  Measurements  as  a Function 
Maximum  Wave  Height  (f  /f  <1) 


Figure  5.5  Fractional  Depth  Loss  Predictions  and  Experimental  Results  as  a 
Function  of  Nonditnensional  Frequency  for  Various  Values  of 
Nondimensional  Input  Displacement  Amplitude 


(a/d)  = 0.4 


5.6  Fractional  Depth  Loss  Predictions  and  Experimental  R 
Function  of  Nondimensional  Frequency  for  Various  Values 
Nondimensional  Input  Displacement  Amplitude 


following  calculations  will  only  consider  fluid  loss  caused  by  horizontal 
ground  motion  and  overpressure  loading.  We  shall  show  that  svch  a predic- 
tion procedure  is  consistent  with  the  field  test  observations  on  MIXED 
COMPANY  [9];  thus  demonstrating  some  confidence  in  applying  this  procedure 
to  estimate  water  loss  In  prototype  situations. 

The  horizontal  low-frequency  ground  shock  environment  on  MIXED 
COMPANY  can  be  approximated  by  a sinusoidal  displacement  waveform  of 
about  one  period  in  duration.  The  basic  frequency  of  the  waveform  was 
about  5 Hz,  and  the  measured  attenuation  of  peak  displacement  with  range 
from  ground  zero  is  given  in  Figure  5.7.  Consequently,  we  shall  model 
D(t)  as  a one-cycle  5-Hz  sinusoidal  waveform  exhibiting  the  attenuation 
charactersitics  given  by  the  line  in  Figure  5.7.  Since  0»5,  the 

predictions  are  relatively  insensitive  to  errors  in  this  nominal  frequency 
value . 

For  the  following  calculations,  the  overpressure  time  history  is 
modeled  as  a sum  of  three  exponentials: 

4P  - P(l-,)(ae-“’  + be'®’  + ce"’'')'' 

4.  4. 

where  t = (t-t  )/D  , t is  the  time  of  shock  arrival  and  D is  the  positive 
s s 

phase  duration.  The  parameters  a,  b,  c,  a,  6,  and  y are  given  as  a func- 
tion of  peak  overpressure  in  Reference  3.  Although  this  overpressure 
approximation  is  valid  only  for  nuclear  environments,  it  can  be  a reason- 
able model  for  HE  explosions  provided  the  actual  HE  overpressures  and 
positive  phase  durations  are  used.  Accordingly,  these  data  are  taken 
from  Reference  11.  For  simplicity,  the  ground  shock  and  airblast  arrival 


Here,  P(x)  and  t^  = x/u,  where  u is  the  approximately  constant 
shock  speed  across  the  containers,  are  the  only  parameters  treated  as 
functions  of  position  over  the  fluid's  surface.  All  other  parameters 
are  taken  to  be  independent  of  position,  bning  functions  of  the  refer- 
ence overpressure  P(x-O),  alone.  This  approximation  is  conservative 
because  it  tends  to  overestimate  the  overpressure  positive  phase  Irjpulse 
gradient  across  the  container.  More  accurate  overpressure  models  are 
currently  being  formulated. 
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times  are  assumed  to  coincide.  Only  a small  error  (usually  <10  percent 
in  Ad/d)  is  introduced  by  this  assumption. 

Figure  5.8  compares  the  small  container  data  with  calculated  water 
loss  from  combined  airblast  and  horizontal  ground  motion  loading.  Vertical 
ground  motion,  debris  impact,  and  dynamic  pressure  effects  are  ignored  in 
this  model.  The  theory  does  differentiate  between  broadside  and  end-on 
configurations;  however,  the  differences  are  too  small  to  indicate  in 
Figure  5.8.  The  broadside  configuration  loses  more  fluid  until  the 
range  exceeds  about  700  ft,  after  which  the  end-on  configuration  loses 
more.  The  calculations  appear  to  qualitatively  describe  the  variation  of 
fluid  loss  with  range,  although  they  generally  underestimate  the  observed 
fluid  loss. 

Comparison  of  calculated  results  and  laboratory  experiments 
with  the  field  experiments  suggests  that  about  half  of  the  observed 
fluid  loss  from  the  small  containers  on  MIXED  COMPANY  can  be  attributed 
to  horizontal  ground  motion  alone , provided  the  range  from  ground  zero  is 
less  than  about  450  ft.  At  ranges  in  excess  of  about  600  ft,  the  theory 
predicts  that  the  airblast  loading  causes  most  of  the  fluid  loss. 


300  400  500  600  700  800  900 

RANGE  - ft 

Figure  5.8  MIXED  COMPANY  Data  and  Predictions  (Neglecting  Vertical  Grdund  Motion,  Debris 
Impact  and  Dynamic  Pressure  Effects)  for  Small  Containers 


The  low-frequency  ground  motions  from  a near-surface  explosion  may 
be  represented  by  a single-cycle  displacement  waveform  whose  amplitude 
is  derived  from 

a^l  0.A5  for  ^ 7.5  (6.1) 

where  V is  the  apparent  crater  volume,  and  R is  the  range  from  ground  zero. 
This  assumption  is  based  on  the  correlation  of  ground  motion  data  from 
existing  high-explosive  and  nuclear  cratering  explosions,  and  is  referred 
to  as  the  "crater  volume"  correlation  within  the  U.S.  ground  shock  com- 
munity. Although  the  data  analysis  leading  to  Equation  6.1  has  not  been 
published  in  the  unclassified  literature,  the  crater  volume  correlation 

has  been  previously  referred  to  in  studies  of  ground  shock  in  hard  rock 
1/3  > 

[12],  For  R/V  'V.  7.5,  Equation  6.1  is  expected  to  underestimate  the  peak 

displacements  which  attenuate  less  rapidly  at  greater  ranges.  This  fact, 

along  with  the  observation  that  ground  motions  become  more  oscillatory 

with  increasing  range,  leads  us  to  believe  that  the  following  results 

1/3  ^ 

underestimate  water  loss  if  R/V  '^7.5. 

8 3 

We  estimate  that  V n.  2 x 10  W ft  for  wet  soils,  were  W is  the  yield 

7 3 

in  megatons.  For  dry  soils,  we  estimate  that  V 5 x 10  W ft  . Since  for 
P ^ 100  psi, 

P 2l  3750 

where  R is  the  range  in  ft.  Equations  6.1  and  6.2  can  be  combined  with  the 
assumption  of  a single-cycle  sinusoidal  waveform  to  obtain 

D (t)  8.3  W^^^(P/600)  sin2iTft  wet  soil,  (6.3) 

w — 

and 

D, (t)  ^1.3  W^^^(P/600)  sin2iift  dry  soil.  (6.4) 


Based  on  analysis  of  ground  motion  data  from  high-explosive  and  nuclear 
cratering  tests  In  the  Pacific,  we  estimajte  that 

f ^ Hz, 

where  0.25 K 1. 

For  the  following  calculations,  the  overpressure,  as  in  the  previous 
chapter,  is  modeled  after  Erode  [3],  and  the  ground  shock  and  air  shock 
arrival  times  are  again  assumed  to  coincide.  Debris  Impact,  vertical 
motion  and  dynamic  pressure  are  completely  ignored. 

6 . 2 Prediction  Calculations 

Figures  6.1  through  6.4  give  predictions  of  fluid  loss  based  on 
Equation  5.4,  where  the  maximum  wave  height  is  calculated  from 
Equation  3.11,  for  nuclear  environments  and  rectangular  contained 
geometries  of  interest.  Figures  6.1  and  6.2  indicate  the  depth  loss 
experienced  under  dry  geologic  conditions  as  a function  of  peak 
overpressure  and  yield.  The  container  dimensions  for  the  calculations 
are  £ = 100  ft,  d = 40  ft  in  Figure  6.1,  and  £ = 300  ft,  d = 40  ft  in 
Figure  6.2.  The  same  Information  for  wet  geologic  conditions  is  presented 
in  Figures  6.3  and  6.4.  The  bands  in  these  figures  arise  from  the 
uncertainty  in  predicting  the  displacement  waveform  frequency.  Generally 
speaking,  the  larger  losses  are  associated  with  the  lower  frequency  limit, 
f = 0.25/W^^^;  the  lower  losses  with  the  higher  frequency  limit, 
f = 1.0/W^^^.  An  exception  to  this  rule  is  the  100  MT,  £ = 100  ft, 
predictions.  In  these  cases,  the  lower  frequency  limit  corresponds  to 
f^/f^l,  while  for  higher  frequency  limit,  f^/f£).  Thus,  an  intermediate 
frequency  actually  is  associated  with  the  larger  loss. 

Because  larger  containers  are  less  sensitive  to  ground  motions  than 
smaller  containers,  the  £ = 300  ft  container  generally  loses  less  fluid 
than  the  £ = 100  ft  container.  Figure  6.2  indicates  that  the  sensitivity 
of  fluid  loss  to  yield  is  extremely  low  for  the  £ = 300  ft  container  at 
low  yields.  In  these  cases,  ground  motion  is  relatively  unimportant,  and 
the  overpressure  is  decaying  quite  rapidly  relative  to  the  normal  mode 
period.  Hence,  Equation  3.13  should  apply,  and  as  mentioned  in  Section  2, 
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Figure  6.1  Predicted  Depth  Loss  as  a Function  of  Peak  Overpressure 
and  Yield  for  Dry  Geologic  Conditions 


gure  6.3.  Predicted  Depth  Loss  as  a Function  of  Peak  Overpressure 
and  Yield  for  Wet  Geologic  Conditions 


the  wave  height  so  determined  ia  effectively  yield-independent.  F)uld 
loss  from  even  large  pools  in  wet  geologies  is  usually  dominated  by  • 

ground  motion. 

The  relative  portion  of  fluid  loss  associated  with  ground  motion 
varies  considerably  with  the  yield  and  container  dimension  variations 
presented  in  Figures  6.1  through  6.4.  For  the  1 “ 300  ft  configuration 
dry  geologic  conditions  and  low  yields,  the  ground  motion  contribution  is 

negligible,  while  the  overpressure  contribution  is  negligible  for  the 

1 - 100  ft  configuration  and  wet  geologic  conditions.  The  comparative 
danger  associated  with  siting  prototype  systems  in  wet  geologies  is 
obvious  from  a comparison  of  Figures  6.1  and  6.3. 

The  fluid  loss  results  presented  in  Figures  6.1  through  6.4,  although 

valid  only  tor  rectangular  containers,  may  be  interpreted  as  a bounding 
I 

estimate  of  fluid  loss  from  any  sloped-wall  container  with  a depth  and 
maximum  horizontal  dimension  given  by  the  same  d and  1 Indicated  in 
these  figures.  This  fact  follows  from  the  observed  wall  slope  dependence 
of  r,  discussed  In  Subsection  5.2. 

Water  loss  from  prototype  trenches,  canals,  etc.,  may  be  prevented, 
or  significantly  reduced  with  the  use  of  berms  placed  around  their  peri- 
meters. As  a practical  matter,  the  cost  of  these  berms  would  be  nominal 
because  excavated  soil  would  be  readily  available  for  their  construction. 
The  wave  height  calculations  discussed  In  this  paper  could  be  used  to 
develop  design  criteria  for  berm  features  (e.g.,  berm  height)  once  the 
threat  environment  and  soil  geology  are  specified. 
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' Provided  vertical  motion,  dynamic  pressure  and  debris  impact  are 
neglected. 


SECTION  7 

CONCLUSIONS  AND  SUMMARY 

7.1  Scaling  Requirements 

Section  2 contains  a discussion  of  scaling  requirements  that  are 
applicable  to  the  design  of  scaled  experiments.  In  the  absence  of  fluid 
compressibility,  viscosity,  surface  tension  and  cavitation,  the  scaling 
laws  for  ground  motion  and  alrblast  overpressure  effects  are  given  by 
Equations  2,7  and  2.8,  respectively.  

Application  of  these  laws  to  high  explosive  field  test  experiments 
presents  difficulties,  however,  since  complete  control  over  the  environ- 
ment is  not  possible.  In  fact,  not  only  are  the  requirements  of 
Equations  2.7  and  2.8  often  mutually  exclusive,  but  la  addition. 

Equation  2.8  cannot  be  satisfied  precisely.  This  observation  led  us 
to  develop  two  approximate  scaling  requirements,  one  valid  in  the  limit 
of  large  containers  and/or  sma7,l  yields  (Overpressure  Impulse  Scaling) 
and  the  other  valid  in  the  limit  of  small  containers  and/or  large  yields 
(Quarter  Root  of  Yield  Scaling).  Impulse  scaling  is  of  more  practical 
interest  since  most  experimental  parameters  fall  into  its  regime  of 
validity. 

7.2  Analytical  Results 

All  analytical  results  presented  here  are  obtained  by  linearizing 
the  appropriate  equations  of  motion  and  boundary  conditions.  This  process 
is  detailed  in  Equations  (3.1)  through  (3.7).  A modal  description  of 
fluid  response  to  horizontal  ground  motion  and  overpressure  loading 
culminates  in  Equation  (3.11).  In  Subsection  3.2  a fast  decay  approxi- 
mation for  the  overpressure  contribution,  (consistent  with  Impulse  Scaling) 
is  derived,  while  in  Subsection  3,3  we  present  a stepwave  approximation 
(consistent  with  Quarter  Root  of  Yield  Scaling) . 

In  Section  4,  we  examine  the  invlscld  coupling  (Kelvin-Helmholtz) 
between  a subsonic  airflow  and  the  fluid.  An  approximate  correction  to 
previously  derived  results  to  account  for  this  interaction  is  contained 


50 


In  Equation  4.7.  It  Is  demonstrated  that  the  magnitude  of  the 
correction  increases  for  higher  mode  numbers  and  large  yield  environ- 
ments. Relatively  speaking,  its  importance  diminishes  as  container 
dimensions  increase.  Also  presented  in  Section  4 is  ah  appfcxlroate 
treatment  of  vertical  ground  motion  effects. 

7.3  Experimental  Results 

A discussion  of  laboratory  ground  shock  simulation  experiments 

is  presented  in  Section  3.  These  experiments  were  conducted  on  a dual 
axis  hydraulic  dynamic  simulator,  and  consisted  of  parametric  studies 
that  varied  input  waveform  frequency,  amplitude,  duration  and  vertical- 
horizontal  phasing.  We  observed  a strong  correlation  between  measured 
depth  loss  and  maximum,  calculated  first  mode  wavehelght,  as  indicated 
in  Figure  5,3.  The  observed  constant  of  proportionality  for  rectangular 
containers  is  approximately  0.3.  Thus  the  data  seem  to  be  well-described 
by 

Ad2i0.3h,  (7.1) 

where  Ad  is  the  depth  loss  and  h is  the  maximum  first  mode  wavehelght. 

Equation  (7.1)  was  assumed  to  hold  for  overpressure  as  well  as 
ground  motion  environments,  and  predictions  for  the  MIXED  COMPANY  fluid 
slosh  experiments  [9]  were  then  generated.  Figure  5.8  compares  the 
MIXED  COMPANY  data  to  the  predicted  depth  losses  given  by  Equation  (7.1). 
Generally  good  agreement  is  observed  although  the  predictions  lie  on  the 
low  side  of  the  data  scatter.  The  same  technique  is  then  used  in 
Section  7 to  estimate  fluid  loss  from  large  scale  containers  exposed  to 
nuclear  environments.  Our  findings  and  conclusions  are  presented  in 
Subsection  6.2. 

7.4  Recommendations 

Technical  areas  where  further  work  is  necessary  for  a complete 
description  of  blast  and  shock  Induced  fluid  motion  include: 

a)  compressible  fluid  effects 

b)  pressure  coupling  between  air  and  fluid 
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shear  Interaction  between  air  and  fluid 
d)  effect  of  debris  impact 

DM  is  currently  studying  the  effect  of  fluid  compressibility.  Pre- 
liininary  Indications  are  that  compressibility  plays  an  inslgnlfleant 
role  in  the  late-tlme  fluid  notion  and  that,  consequently,  the  results 
described  above  will  not  be  significantly  affected.  No  such  tentative 
conclusion  can  be  made,  however,  for  the  remaining  effects  listed  above 
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LIST  OF  SYMBOLS 


peak  displacement  amplitude  of  horizontal  ground  motion 

modal  coefficients  for  velocity  potential 

modal  coefficients  for  surface  displacement 

drag  coefficient 

undisturbed  fluid  depth 

horizontal  ground  motion  displacement 

vector  ground  motion  displacement 


- positive  phase  duration 

- temporal  frequency 

- frequency  structure  function  for  maximum  wave  height 

- gravitational  constant 

- maximum  first  mode  wave  height 

- vertical  ground  motion  displacement 

- unit  vector  in  x-direction 

- positive  phase  overpressure  impulse 

- surface  disturbance  wave  number 

- unit  vector  in  z direction 

- container  length 

- mach  number  of  airflow 

- mode  number 

- alrblast  overpressure 

- P/pg£ 

- alrblast  dynamic  pressure 

- radius  from  ground  zero 


slosh  period 
i/u 


- air  shock  propagation  velocity 

- crater  volume 
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LIST  OF  SYMBOLS  (continued) 


- explosive  yield 

- horizontal  dimension 

- vertical  dimension 


modal  coefficients  for  ground  displacement 

modal  coefficients  for  alrblast  overpressure 

overpressure  perturbation  induced  by  free  surface  disturbance. 

fluid  depth  loss 

radian  frequency  of  n^^  mode 

time  decay  constant  for  dynamic  pressure 

fluid  velocity  potential 

reduced  velocity  potential 
fluid  density 

constant  of  proportionality  between  wave  height  and  depth  loss 

t//T/g~ 

gradient  operator 


free  surface  displacement 

n/J. 
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ATTN:  John  Collins 


J.  L.  Merritt 

Consulting  & Special  Engr,  Svs, , luo, 
ATTN:  J.  L.  Merritt 
ATTN:  Tech.  Lib. 


Kaman  AviDyne 

Division  of  Kaman  Sciences  Corporation 


ATTN 

ATTN 

ATTN 

ATTN 


Norman  P.  Hobbs 
E,  S.  Crisolone 
O.  Zartarlan 
Tech.  Lib. 


Kaman  Sciences  Corporation 


ATl'N 

ATTN 

ATTN 


Paul  A.  Ellts 
Frank  H.  Shelton 
Library 


Lockheed  Missiles  ti  Space  Co.  Inc, 
ATTN;  Tech.  Lib. 


Lockheed  Missiles  A Space  Company 

ATTN;  Tom  Geers,  Dept.  62-33,  Bldg.  205 
ATTN;  Tech.  Info.  Ctr. , D/Coll. 


Lovelace  Foundation  for  Medical  Education 

ATTN;  Aast.  Dtr.  of  Res. , Robert  K,  Jones 
ATTN;  Tech.  Lib. 


Martin  Marietta  Aerospace 
Orlando  Division 

ATTN;  G,  Fotteo 


University  of  Maryland 
Dept,  of  Civil  Engineering 
ATTN:  Bruce  S.  Berger 


McDonnell  Douglas  Corporation 
ATTN:  Robert  W,  Halprln 


Meteorology  Research,  Inc. 
ATTN:  William  D,  Green 


The  Mitre  Corporation 
ATTN;  Library 


Nathan  M.  Newmnrk 
Consulting  Engineering  Services 
ATTN;  Nathan  M,  Newmnrk 


PacUtca  Technology 
ATTN;  J.  Kent 
ATTN;  H.  BJork 


Physics  International  Company 


ATTN 

ATTN: 

ATTN 

ATTN: 

ATTN 

ATTN 

ATTN 

ATTN 


Doc,  Con.  for  Dennis  Orphal 
Doc.  Con.  for  E,  T.  Moore 
Doc.  Con,  tor  Larry  A.  Behrmann 
Doc.  Con.  for  Tech.  Lib. 

Doc.  Con.  tor  Coyo  Vincent 
Doc.  Con.  for  Fred  M.  Sauer 
Doc.  Con.  for  Robert  Swift 
Doc.  Con.  for  Charles  Godfrey 
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FolytechnlQ  Inst,  of  Brooklyn 
Dept,  of  Aerospace  & Applied  Mach. 

ATTN!  J.  M.  Kloaner 

R & D Associates 

ATTN:  William  B,  Wright,  Jr. 

ATTN:  Sheldon  Schuster 
ATTN:  Tech,  Lib. 

ATTNi  Albert  L,  Latter  

ATTNt  Cyrus  P.  Knowles 
ATTN:  Bruce  Hartenbaum 
ATTN:  J.  G.  Lewis 
ATTN:  Henry  Cooper 
ATTN:  Harold  L.  Erode 
ATTN:  Jeriy  Carpenter 
ATTN:  E.  J.  Chapyak 

The  Rand  Corporation 
ATTN:  C.  C.  Mow 

Research  Analysis  Corporation 
ATTN:  Documents  Library 

Science  Applications,  Inc. 

ATTN:  D.  E.  Maxwell 
ATTN:  David  Bernstein 

Science  Applications,  Inc. 

ATTN:  William  M,  Layson 
ATTN:  R.  Seebaugh 
ATTN:  John  Mansfield 

Science  Applications,  Inc. 

ATTN;  Michael  McKay 
ATTN;  Tech.  Lib. 

Science  Applications,  Inc, 

ATTN:  R.  A,  Shunk 

Shock  Hydrodynamics,  Inc. 

A Dlv.  of  Whittaker  Corporation 
ATTN:  L,  Zemow 

Southwest  Research  Institute 
ATTN:  Wilfred  E.  Baker 
ATTN:  A,  B.  Wenzel 

Stanford  Research  Institute 

ATTN:  SR[  Library,  Rm,  G-021 
ATTN;  Carl  Peterson 
ATTN:  BurtR.  Gasten 
ATTN:  George  R.  Abrahamson 

Sundstrand  Corporation 

ATTN:  Curtis  B.  White 


DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 

Systems,  Science  It  Software,  Ino. 

:ATTN:  Ted  Cherry  1 T 

ATTNi  Robert  T.  Allan 
ATTN:  Thomas  D.  Blnoy 
ATTNi  Donald  R.  Orine 
ATTNi  Tech.  Lib. 

Terra  Tek,  Inc. 

ATTNi  A.  H.  Jones 
" ATTNi  Sidney  Green 
ATTNi  Tech.  Lib. 

Tetra  Tech,  Ino, 

ATTNi  LI -San  Hwang 
ATTNi  Toch.  I-b. 

TRW  Systems  Grou,. 

ATTNi  Benjamin  Suasholtz 

ATTNi  Norm  LIpner 

ATTNi  William  Rowan 

ATTNi  Jack  Farrell 

ATTNi  Tech.  info.  Ctr. , S-I930 

ATTNi  Pravln  Bhutta 

ATTNi  Paul  Lteberman 

TRW  Systems  Group 
San  Bernardino  Operattona 
ATTNi  GregHulcher 

Universal  Analytics,  Inc. 

ATTNi  E.  1.  Field 

URS  Research  Company 
ATTNi  Ruth  Schneider 
ATTNi  Tech.  Lib. 

The  Eric  H.  Wang 

Civil  Engineering  Research  Facility 
ATTNi  Larry  Btckle 
ATTNi  Neal  Baum 

Washington  State  Unlversl^ 

Administrative  Office 

ATTNi  Artltur  Miles  Hohoif  for  George  Duval 

Weidllnger  Assoc.  Consulting  Engineers 
ATTNi  Melvin  L.  Baron 

Weidllnger  Assoc.  Consulting  Engineers 
ATTNi  J.  Isenberg 

V'estlnghouse  Electric  Company 
Marine  Division 

ATTNi  W.  A.  Volz 


«2 


